Constraints on Scalar Spectral Index from Latest Observational Measurements 
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Recently, the nine-year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9) 
found that the inflationary models with the scalar spectral index > 1 are excluded at about 
5 a confidence level. In this paper, we set the new limits on the scalar spectral index in different 
cosmological models combining the WMAP9 data with the small-scale cosmic microwave background 
measurement from the South Pole Telescope, baryon acoustic oscillation data, Hubble Telescope 
measurements of the Hubble constant, and supernovae luminosity distance data. In most of extended 
cosmological models, e.g. with a dark energy equation of state, the constraints on Us do not change 
significantly. The Harrison-Zel'dovich-Peebles (HZ) scale invariant spectrum is still disfavored at 
more than 4 a confidence level. However, when considering the model with a number of relativistic 
species A^off , we obtain the limit on the spectral index of ris — 0.980 ± 0.011 (1 a), due to the strong 
degeneracy between Hs and Ngff. The HZ spectrum now is still consistent with the current data at 
95% confidence level. 



I. INTRODUCTION 

Inflation, the most attractive paradigm in the very 
early universe, has successfully resolved many problems 
existing in the hot big bang cosmology, such as the 
flatness, horizon, monopole problem, and so forth [1]. 
Its quantum fluctuations turn out to be the primor- 
dial density fluctuations vifhich seed the observed large 
scale structures (LSS) and the anisotropics of cosmic mi- 
crowave background (CMB). To distinguish various infla- 
tionary models, the spectral index of the power spectrum 
of primordial curvature perturbations is one of the most 
important variables. 

With the accumulation of observational data from 
CMB, LSS and Type la Supernovae observations (SN) 
and the improvements of the data quality, the cosmolog- 
ical observations play a crucial role in our understanding 
of the Universe and also in constraining the cosmological 
parameters [2-5]. Thus, determining the scalar spectral 
index from the observational data is a very powerful 
and reliable way to understand these inflationary models. 

Recently, the new CMB data have been released 
[6, 7]. The nine-year data release of Wilkinson Mi- 
crowave Anisotropy Probe (WMAP9) has determined 
the cosmological parameters accurately and found that 
the 68% C.L. constraint on the scalar spectral index of 
= 0.9608 ±0.0080 [7], when combining with the small- 
scale CMB measurement from the South Pole Telescope, 
baryon acoustic oscillation data (BAO) and Hubble Tele- 
scope measurements of the Hubble constant (HST). 
Within the ACDM framework, the Harrison-Zel'dovich- 
Pcebles (HZ) scale invariant spectrum (??s = 1) and the 
spectra with > 1 arc disfavored by about 5 a confi- 
dence level. 

Although the ACDM model is a good candidate for 
interpreting the data, the evidence for various extended 



models should not be neglected, such as the number of 
relativistic species [8, 9], the massive neutrino [10, 11] 
or the equation of state of dark energy [12]. More im- 
portantly, the degeneracies between ris and cosmological 
parameters introduced in these models could weaken the 
constraint on ng [13-15]. 

In this paper, we explore the cosmological constraints 
on Us in some extended models from the latest cosmolog- 
ical data sets, including the WMAP9 temperature and 
polarization power spectra, the small-scale CMB mea- 
surement from SPT, the BAO measurements from sev- 
eral LSS surveys, the HST prior on the Hubble constant 
Hq and the "Union2.1" compilation SN sample made by 
the Supernova Cosmology Project. Firstly, wc consider 
the general inflationary model with the tensor fluctua- 
tions (r) in the ACDM framework. We then extend the 
ACDM model allowing for the dark energy models with 
a constant equation of state (EoS, w) or with a time- 
varying EoS {w{z)). Finally, we include the massive neu- 
trino case (Y] ^1^) or the number of relativistic species 
(iVeff) in the ACDM model. 

Our paper is organized as follows: In Section II we de- 
scribe the method and the latest observational data sets 
used in the numerical analyses; Section III contains our 
main global constraints of the scalar spectral index Us in 
different cosmological models from the current observa- 
tions. The last Section IV is the conclusions. 



II. METHOD AND DATA 

A. Numerical Method 

We perform a global fitting of cosmological parame- 
ters using the CosmoMC package [16], a Markov Chain 
Monte Carlo (MCMC) code. We assume purely adi- 
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abatic initial conditions and a flat ACDM Universe. 
The following six cosmological parameters are allowed 
to vary with top-hat priors: the cold dark matter en- 
ergy density parameter flch'^ e [0.01,0.99], the baryon 
energy density parameter n^h^ £ [0.005,0.1], the scalar 
spectral index Us G [0.5,1.5], the primordial amplitude 
ln[10i°A^] € [2.7,4.0], the ratio (multiplied by 100) of 
the sound horizon at decoupling to the angular diameter 
distance to the last scattering surface 1009s G [0.5, 10], 
and the optical depth to reionization t e [0.01, 0.8]. The 
pivot scale is set at fc^o = 0.05Mpc~^. Besides these 
six basic cosmological parameters, we have several ex- 
tra cosmological parameters in different extended cos- 
mological models: the running of scalar spectral index 
as = dlnus/dlnk £ [—0.1,0.1]; the tensor to scalar ra- 
tio of the primordial spectrum r = At/Ag G [0,2]; the 
total neutrino mass fraction at the present day 



VLmh? ~ 93.14 eV VLmh? 



e [0,0.1] 



(1) 



and the number of relativistic species iVcff e [0, 10]. We 
also consider the dark energy model with the EoS param- 
eters Wo £ [—2, 0] and wi £ [—5, 2], which is given by the 
parametrization [17] 



Wrfe(a) = Wo + wi(l - a) 



(2) 



where a = 1/(1 -|- z) is the scale factor and wi ~ —dw/da 
characterizes the "running" of EoS. The ACDM model 
has Wo = —1 and wi = 0. For the dark energy model 
with a constant EoS, wi — 0. When using the global fit- 
ting strategy to constrain the cosmological parameters, 
it is crucial to include dark energy perturbations [18]. In 
this paper we use the method provided in refs. [18, 19] 
to treat the dark energy perturbations consistently in 
the whole parameter space in the numerical calculations. 
Therefore, the most general parameter space in the anal- 
yses is: 

{nbh'^,Qch'^,Qs,T,ns,As,as,r,WQ,wi,f^,Ncs} ■ (3) 



B. Current Observational Data 

In our analysis, we consider the following cosmological 
probes: i) power spectra of CMB temperature and po- 
larization anisotropies; ii) the baryon acoustic oscillation 
in the galaxy power spectra; iii) measurement of the cur- 
rent Hubble constant; iv) luminosity distances of type la 
supernovac. 

To incorporate the WMAP9 CMB temperature and 
polarization power spectra, we use the routines for com- 
puting the likelihood supplied by the WMAP team 
[7]. The WMAP9 polarization data are composed of 
TE/EE/BB power spectra on large scales (2 < £ < 23) 
and TE power spectra on small scales (24 < £ < 800), 
while the WMAP9 temperature data includes the CMB 
anisotropics on scales 2 < £ < 1200. Furthermore, we 



also use the recent SPT data [20], using 47 bandpowcrs 
in the range 600 < £ < 3000. The likelihood is assumed 
to be Gaussian, and we use the published band-power 
window functions and covariance matrix. In order to 
address for foreground contributions, the SZ amplitude, 
the amplitude of the clustered point source contribution, 
and the amplitude of the Poisson distributed point source 
contribution, are added as nuisance parameters in the 
CMB data analyses. 

Baryon Acoustic Oscillations provides an efficient 
method for measuring the expansion history by using 
features in the clustering of galaxies within large scale 
surveys as a ruler with which to measure the distance- 
redshift relation. It provides a particularly robust quan- 
tity to measure [21]. It measures not only the angular 
diameter distance, Da{z), but also the expansion rate of 
the universe, H(z), which is powerful for studying dark 
energy [22]. Since the current BAO data are not accu- 
rate enough for extracting the information of DaIz) and 
H{z) separately [23], one can only determine an effective 
distance [24]: 



D,iz) = [il + zfD\iz)cz/Hiz)]'/^ 



(4) 



In this paper we use the recent BAO measurement at high 
rcdshift z = 2.3 detected in the Ly-a forest of Baryon Os- 
cillation Spectroscopic Survey (BOSS) quasars [25]. Fur- 
thermore, we also include the BAO measurement from 
the 6dF Galaxy Rcdshift Survey (6dFGRS) at a low red- 
shift z ~ 0.106 [26], and the BAO measurements from the 
WiggleZ Survey at three rcdshift bins z — 0.44, z = 0.60 
and z = 0.73 [27], the measurement of the BAO scale 
based on a re-analysis of the Luminous Red Galaxies 
(LRG) sample from Sloan Digital Sky Survey (SDSS) 
Data Release 7 at the median rcdshift z ~ 0.35 [28], 
and the BAO signal from BOSS CMASS DR9 data at 
z = 0.57 [29]. 

In our analysis, we add a Gaussian prior on the cur- 
rent Hubble constant given by ref. [30]; Hq = 73.8 ± 2.4 
kms^^Mpc^^ (68% C.L.). The quoted error includes 
both statistical and systematic errors. This measurement 
of Hq is obtained from the magnitude-redshift relation of 
240 low-z Type la supernovac at z < 0.1 by the Near 
Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) Camera 2 of the Hubble Space Telescope (HST). 
This is a significant improvement over the previous prior, 
i/o = 72 ± 8 kms"^Mpc"\ which is from the Hubble 
Key project final result. In addition, we impose a weak 
top-hat prior on the Hubble parameter: Hq £ [40, 100] 
kms~^ Mpc"-'^. 

Finally, we include data from Type la supernovac, 
which consists of luminosity distance measurements as a 
function of rcdshift, Dl{z). In this paper we use the lat- 
est SN data sets from the Supernova Cosmology Project, 
"Union Compilation 2.1", which consists of 580 samples 
and spans the rcdshift range < z < 1.55 [31]. This 
data set also provides the covariance matrix of data with 
and without systematic errors. In order to be conserva- 
tive, we use the covariance matrix with systematic errors. 
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FIG. 1: Marginalized one-dimensional and two-dimensional 
likelihood (1, 2 a contours) constraints on the parameters Us 
and Ho in the standard ACDM model hom different present 
data combinations: WMAP9 only (red), WMAP9-I-SPT 
(blue) and All datasets (green). 



When calculating the likelihood from SN, we marginal- 
ize over the absolute magnitude M, which is a nuisance 
parameter, as done in refs. [32]. 

III. NUMERICAL RESULTS 

In this section we present our global fitting results of 
the cosmological parameters determined from the latest 
observational data and focus on the degeneracies between 
Us and other extended parameters in different models. 

A. Standard ACDM model 

Firstly, we consider the standard ACDM model. In ta- 
ble I we show the constraints on some related cosmologi- 
cal parameters from three different data combinations: 



TABLE I: 1 (7 constraints on some cosmological parameters 
from different data combinations in the standard ACDM 
model. 



Standard ACDM 





WMAP9 alone 


WMAP9+SPT 


All Datasets 




0.972 ±0.013 


0.966 ±0.011 


0.963 ± 0.008 


Ho 


70.34 ±2.21 


70.20 ± 2.15 


68.90 ± 0.62 




2.270 ±0.050 


2.234 ±0.042 


2.224 ± 0.034 




11.37 ±48 


11.41 ±0.47 


11.70 ±0.17 



WMAP9 alone, WMAP9+SPT, and All datasets. In 
the upper panel of figure 1 we show the one-dimensional 
marginalized likelihood distributions of ??s from three 
data combinations. Using the WMAP9 data alone, we 
obtain the 68% constraint of = 0.972 ±0.013. The pri- 
mordial spectra with > 1 are only excluded at 2 cr con- 
fidence level. When we include the small-scale SPT mea- 
surement, the constraint on becomes slightly tighter. 
Us = 0.966 ± 0.011 at 1 cr confidence level. Since the 
median value and the error bar of rig are smaller, when 
comparing with those from WMAP9 alone, the signifi- 
cance of ris < 1 is more than 3 a confidence level. We 
also show the two-dimensional contour between Hq and 
ns in the below panel of figure 1. There is a strong cor- 
relation between and Hq, when considering the CMB 
data alone. 

When combining all the datasets together, this strong 
degeneracy is broken apparently. The constraint on the 
Hubble constant becomes much more stringent, ho = 
0.6890 ± 0.0062 (Icr C.L.). Consequently, the constraint 
of the spectral index also becomes tighter significantly. 
Us ^ 0.963 ± 0.008 (Icr C.L.). The error bar of is 
reduced by a factor of 1.5, due to the constraining power 
of BAG, HST and SN. The HZ spectrum is disfavored by 
the current data at about 5cr confidence level, which is 
consistent with that from the WMAP9 paper [7]. 

However, this strong constraint on the spectral index is 
model-dependent. In some extended ACDM models, the 
constraints on Ug could be changed, due to the possible 
degeneracies between and other extended parameters 
in different models. In the following sections, we discuss 
the constraints on parameters in these extended cosmo- 
logical models, which is shown in table II, such as the 
inflationary models, as and r, the dynamical dark energy 
model, w{z), and ones including the neutrino properties, 
X] and Neff. 

B. Inflationary Models 

Firstly, we include the gravitational waves into the 
analysis. When using all datasets together, the data yield 
the 95% upper limit of tensor-to-scalar ratio r < 0.15. 
Meanwhile, the constraint on the spectral index is slightly 
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FIG. 2: Marginalized two-dimensional likelihood (1, 2(7 con- 
tours) constraints on the parameters ris and r from all 
datasets together (red). The two blue solid lines are predicted 
by the m?4>^ and A(^* models, respectively. The green points 
denote predictions assuming that the number of e-foldings 
Af = 50 — 60 from two models. 



TABLE II: 1(7 constraints on cosmological parameters r, Qs, 
rriv and dark energy EoS from the current observations in 
different extended models. For the weakly constrained pa- 
rameters we quote the 95% upper limits instead. 



models 


constraints 


Us constraints 


ACDM-fr- 


r < 0.15 


0.966 ± 0.009 


ACDM+Os 


as = -0.023 ±0.011 


0.948 ±0.011 


WCDM 


w = -1.060 ±0.066 


0.960 ± 0.009 


W(z)CDM 


wo = -1.11 ± 0.15 


0.962 ±0.011 




wi =0.18 ±0.65 




ACDM-fX^'"-- 


™^ < 0.47eV 


0.968 ± 0.009 



relaxed, Ug = 0.966 ± 0.009 at 68% confidence level, 
due to the degeneracy between and r. In figure 2 
we show the two-dimensional constraints in the {ris^r) 
panel which can be compared with the prediction of the 
inflation models. We find that the HZ scale-invariant 
spectrum (n^ = 1, r = 0) is still disfavored at about 4cr 
confidence level. Also, the inflation models with "blue" 
tilt {us > 1) are excluded by the current observations. 
Furthermore, assuming the number of e-foldings N = 
50 — 60, the single slow-rolling scalar field with poten- 
tial V{(t)) - m202, which predicts (ns,r)=(l -2/N,8/N), 
is still within the 2 a region, while another single slow- 
rolling scalar field with potential V{(t)) ^ A(^^, which pre- 
dicts {ns,r) = {l — 3/N,16/N), has been excluded more 
than at 2 cr confidence level. 

We also explore the constraint on the running of the 
spectral index from the latest observational data. When 
neglecting the tensor fluctuations (r = 0), the combina- 
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FIG. 3: Marginalized two-dimensional likelihood (1, 2(7 con- 
tours) constraints on the parameters ris and as from all 
datasets together with (blue) and without (red) considering 
the tensor fluctuations in the analysis. 



tion of the current observational data yield the limit on 
the running of the spectral index of as — —0.023 ±0.011 
(Icr), which means the running of Ug is favored by the 
current data at more than 2 a confidence level. In fig- 
ure 3 we show the two-dimensional constraints in the 
(ns,as) panel. Due to the degeneracy between Ug and 
Ug, the 68% constraint on Ug is slightly enlarged, Ug — 
0.948 ±0.011. The error bar is relaxed by a factor of 1.5, 
when comparing with the standard ACDM model. 

Finally, we vary the ag and r simultaneously in the 
analysis. From the blue contour of figure 3, one can 
see that the constraint on Ug does not change, Ug = 
0.949±0.011 (1 a). The degeneracy between ag and r sig- 
nificantly weakens the constraints on them, namely the 
68% constraint on Ug is ag = -0.039±0.016 and the 95% 
upper limit on r is r < 0.35. The current data still favor 
the running of Ug at more than 2 a confidence level. 



C. Dynamical Dark Energy 

Assuming the fiat universe, first we explore the cos- 
mological constraints in the dark energy model with a 
constant EoS, w {w = wq, wi = Q), from the latest ob- 
servational data. In figure 4 we show the two-dimensional 
constraints on w and Ug. Current observational data 
yield a strong constraint on the constant EoS of dark 
energy, w = —1.060 ± 0.066 (Icr), which is similar with 
the limit from WMAP9 [7]. The standard ACDM model 
(w — —1) is consistent with the current observational 
data. In this case the constraint on Ug is slightly changed, 
Ug = 0.960 ± 0.009 at 68% confidence level, due to the 
degeneracy between Ug and w. 

For the time evolving EoS, in figure 5 we illustrate the 
constraints on the dark energy parameters wq and wi. 
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FIG. 4: Marginalized two-dimensional likelihood (1, 2 a con- 
tours) constraints on the parameters ris and w from all 
datasets together (red). 



FIG. 6: Marginalized two-dimensional likelihood (1, 2(7 con- 
tours) constraints on the parameters Us and ^ m,v from all 
datasets together (red). 
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FIG. 5: Marginalized two-dimensional likelihood (1, 2cr con- 
tours) constraints on the parameters wq and wi from all 
datasets together (red). The blue solid lines stand for wq = 
— 1 and Wo -\- — —1. 



FIG. 7: Marginalized two-dimensional likelihood (1, 2a 
contours) constraints on the parameters Us and A'^off from 
different present data combinations: WMAP9 only (red), 
WMAP9+SPT (blue) and All datasets (green). 



For the flat universe, due to the limits of the precisions 
of observational data, the variance of ujq and wi are still 
large, namely, the 68% constraints on wq and wi are 
Wo = -1.11 ± 0.15 and wi = 0.18 ± 0.65. And the 95% 
constraints are —1.38 < wq < —0.80 and —1.32 < wi < 
1.15. This result implies that the dynamical dark energy 
models are not excluded and the current data cannot 
distinguish different dark energy models decisively. The 
obtained best fit model is the Quintom dark energy model 
[35] with the particular feature that its EoS can cross the 
cosmological constant boundary smoothly. The standard 
ACDM model, however, is still a good fit right now. 
For the spectral index, the weak correlations between 



Us and the parameter of dark energy EoS do not change 
the limit significantly, namely the 68% constraint is rig = 
0.962±0.011. The HZ spectrum is stiU ruled out at about 
4 (T confidence level. 



D. Neutrino Properties 

Finally, we study the constraint on Ug in the cosmolog- 
ical models including the neutrino properties, the mas- 
sive neutrino and the number of relativistic species. In 
tabic II wc show the constraint on the total neutrino mass 
from all datasets together, namely the 95% upper limit is 
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TABLE III: 1 a constraints on some cosmological parameters 
from different data combinations in the model with the num- 
ber of relativistic species A'cff . 



ACDM + N^tf 





WMAP9 alone 


WMAP9+SPT 


All Datasets 


Us 


0.994 ± 0.026 


0.986 ±0.018 


0.980 ±0.011 


iVoff 


4.88 ± 2.08 


3.82 ±0.57 


3.79 ±0.38 


Ho 


76.84 ±8.18 


74.31 ± 3.60 


72.29 ± 1.82 


lOOOt/i^ 


2.259 ± 0.050 


2.271 ± 0.048 


2.252 ±0.036 




14.75 ± 3.79 


12.68 ± 1.04 


13.04 ± 0.73 



^ mn < 0.47eV, which is consistent with previous works 
[7, 8, 10, 33, 34]. In figure 6 we show the constraints in 
the {ns,J2''^i^) panel. Since the degeneracy between 
and is not very strong, including the massive neu- 

trino does not change the constraint on rig significantly, 
namely rig = 0.968 ± 0.009 at 68% confidence level. 

Then we consider the constraints on the number of 
relativistic species, A'cff , from different data combi- 
nations (Table III), assuming massless neutrinos. We 
find the WMAP9 data alone gives very weak constraint 
A'cff ~ 4.88 ±2.08 at the 68% confidence level, consistent 
with the result derived by the WMAP9 team. Adding 
the small-scale SPT data significantly improves the con- 
straints on A'^eff to A'^eff = 3.82±0.57 (1 a). When we com- 
bine all datasets together, we obtain our most stringent 
constraint of A^cff = 3.79 ± 0.38 (68% C.L.). Our results 
are quite consistent with previous works [8, 10, 20], and 
display a slight preference for an extra relativistic relic. 
However, the standard value of A'eff = 3.04 remains well 
within the 95% confidence intervals. 

Since A'cff can be written in terms of r2,„/i^ and the red- 
shift of matter-radiation equality, Zcq, there are strong 
degeneracies present between A'cff, the matter density, 
flmh^ and the Hubble parameter Hq. Consequently, 
A'cff is also strong correlated with the spectral index 
rig. In figure 7 we show the constraints on rig and 
A'cff from different data combinations. When using all 
datasets together, the 68% C.L. constraint on rig becomes 
rig = 0.980 ± 0.011. The HZ spectrum now is consistent 
with the current data at 95% confidence level. 

When we include the number of relativistic species, 
A'cff, and the massive neutrino simultaneously, the con- 
straints on parameters become weaker, due to the de- 
generacies among them. The 95% upper limit of the to- 
tal neutrino mass is "Ym^ < 0.71. The 68% C.L. con- 
straints on A'cff and rig becomes A'cff = 3.89 ± 0.39 and 
rig = 0.988 ±0.013. 



IV. SUMMARY 

Recently many experimental groups have published 
their new observational data, such as temperature and 
polarization power spectra of WMAP9 [7], temperature 
power spectrum of SPT at high multipoles £ [20], and the 
BAO measurement from the Ly-a forest of BOSS quasars 
at high redshift z = 2.3 [25]. The WMAP collaboration 
has presented the cosmological implications of their final 
nine- year data release, finding that the spectra with the 
spectral index rig > 1 are disfavored by the current obser- 
vational data at about 5 a confidence level in the ACDM 
framework. 

However, in the analyses we find that the strong con- 
straint on rig could be weakened by considering the pos- 
sible degeneracies between rig and other cosmological pa- 
rameters introduced in some extended models, such as 
the tensor fluctuation r and the dark energy EoS w. The 
largest effect is shown in the model with the number of 
relativistic species A^eff- Due to the strong degeneracy 
between rig and A'cff, the error bar of rig is significantly 
enlarged, namely rig = 0.980 ± 0.011 (la), and the HZ 
spectrum now is consistent with the current data at 95% 
confidence level. 

Finally, we vary all of the extended cosmological pa- 
rameters at the same time. Due to the degeneracies 
among these parameters, their constraints become much 
weaker. More importantly, the limit on the spectral index 
rig is relaxed further, rig — 0.977 ± 0.024 at 68% confi- 
dence level. The disagreement between the HZ spectrum 
and the current observational data disappears. The high 
accurate Planck measurement [36] is needed, in order to 
constrain rig and distinguish various inflationary models 
better. 
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